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We study a waveguide array with an embedded nonlinear saturable impurity. We solve the impurity problem 
in closed form and find the nonlinear localized modes. Next, we consider the scattering of a small-amplitude 
plane wave by a nonlinear impurity mode, and discover regions in parameter space where transmission is fully 
suppressed. We relate these findings with Fano resonances and propose this setup as a mean to control the 
transport of light across the array. © 2009 Optical Society of America 
OCIS codes: 190.0190, 190.5530, 190.6135, 230.4320 



Discrete nonlinear systems have been studied in the 
context of nonlinear optical waveguide arrays during the 
last twenty years [1,2]. One of its main advantages is the 
possibility to change and control all essential parameters, 
e.g. geometry, dimensionality, nonlinearity, beam angle, 
etc. On the other hand, scattering problems have always 
been of great interest to the physics community due to its 
capacity to probe the physical nature of many complex 
systems. In the context of nonlinear waveguide arrays, 
the main properties of soliton scattering against impu- 
rity potentials have been theoretically [3] and experimen- 
tally [4] studied, and a rich and complex phenomenology 
has been observed. During the last years the scattering 
of plane waves by nonlinear localized modes has opened 
the possibility to observe Fano resonances [5] in non- 
linear waveguide arrays: Nonlinearity generates several 
scattering channels which can lead to resonances due to 
destructive interference and, as a consequence, to total 
absence of transmission similar to the original Fano prob- 
lem. As examples we can cite the prediction of Fano res- 
onances in the context of nonlinear quadratic waveguide 
arrays [6] and also, in a very different research area, the 
recent prediction of similar resonances in Bose-Einstein 
condensates [7] (see Ref. [8] as a review of recent findings 
in this area). 

In the present letter we propose a saturable impurity 
embedded in a linear waveguide array as a new possi- 
ble experimental setup to observe Fano resonances. Re- 
cently, scattering of plane waves by bright and dark soli- 
tons has been experimentally studied in saturable non- 
linear media [9] but, to the best of our knowledge, no 
direct observation of this type of resonances has been 
implemented yet. We also characterize the main proper- 
ties of nonlinear impurity modes (NLM) in this type of 
system and suggest the possibility of a switching-mode 
scheme, based on a judicious tuning of the system pa- 
rameters. 

We consider the saturable discrete nonlinear 
Schrodinger (s-DNLS) equation in a dimensionless 
form with a linear defect and nonlinear interaction only 



at one site (n = m) of a ID waveguide array: 
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These coupled-mode equations describe the propagation 
of light in weakly-coupled waveguides close to the first 
band of the band-gap structure. Amplitude u n repre- 
sents the light amplitude at a guide centered on site n, 
z the propagation coordinate along the waveguides, e 
is the strength of the linear defect and (3 the nonlinear 
coefficient. We include both type of defects, linear and 
nonlinear, in order to deal with a more general problem. 
A linear site impurity can be created by altering the ge- 
ometry of a given waveguide, while at the same time 
tuning the spacing with its nearest-neighbors in order to 
keep its coupling with the rest of the array unaltered. On 
the other hand, the nonlinear response at the impurity 
site can be boosted by a judicious amount of extra metal 
doping [10]. 

First, we look for stationary localized solutions cen- 
tered at site rii. We insert into (1) the ansatz u n (z) = 
U n (z) = U x\ n - ni \e iXz with U ,x G R. U is the 
impurity-mode amplitude and x determines its localiza- 
tion length. By defining g = (3/(1 + U§) and imposing 
\x\ < 1 we obtain: x = (g - e)/2 ± y/l + ((g - e)/2) 2 
and A = ±a/4 + (g — e) 2 . The relation between g and e 
determines the sign of A. For g > e (g < e) the sign is mi- 
nus (plus), x and A < (x and A > 0), and the solution 
is staggered (unstaggered). The optical power is defined 
as P = E n K| 2 = U§^4 + (g-e) 2 /\g - e|. Hereafter, 
we fix (3 = 10 and e = 5 (Thus, if U = 1 g = e). Fig- 
ure 1 shows power vs propagation constant of the impu- 
rity mode. From Fig. 1(a) we see that there is no power 
threshold for staggered modes (which are always linearly 
stable). On the other hand, unstaggered modes posses a 
threshold that separates two regimes: dP/dX < un- 
stable and dP/dX > stable solutions. When nonlinear 
modes approach the linear band (\g — e\ — > 0), the optical 
power diverges. This amplitude-dependent "transition" 
implies that, unlike other impurity type, for a single sat- 
urable impurity, one can switch from a staggered to an 
unstaggered mode just by varying the input power. 

Now, we look for Fano resonances, i.e. what would 
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Fig. 1. (a) P vs A diagram for NLM. (b) A vs Uq for 
NLM (black lines) and uo cc vs Uq for local modes (gray 
lines). 

happen, if a small-amplitude plane wave is scattered by a 
nonlinear impurity mode? In the impurity region, where 
the interaction takes place, a "local mode" can be gener- 
ated. This mode corresponds to an extra channel, which 
is opened by the interaction. When this mode is fully 
excited, a zero transmission of plane waves can occur 
due to destructive interference and a Fano resonance ap- 
pears [8]. We assume the plane wave amplitude <j) n (z) 
to be much smaller than the impurity- mode amplitude, 
i.e. 1 0o | <C \Uq\. We insert u n (z) = U n (z) + <j> n {z) an d 
linearize (1) with respect to <j) n (z) obtaining: 



dz 



>n-l + ( e - — ) 0O^n,r 



9 -U 2 e 2iXz ^n, ni . (2) 

To solve this problem, we use the ansatz <j) n {z) = 
a n e luJZ + b n e^ 2X ~^ z in (2), obtaining two coupled dis- 
crete equations: 

ioa n = (a n+ i + a n _i) + La n £ n , n ., (3) 
(2A - u))b n = (6 n+ i + b n -i) + Lb n 5 n , ni , (4) 

where Lc n = (e - g 2 / 0) c n + (g 2 U% / 0) d n with c n 7^ d n 
(c ni d n = a n ,b n ). a n corresponds to the open channel, 
i.e. a traveling plane wave with a frequency given by 
ujk = 2 cos k. k is the plane wave vector which, in a ex- 
periment, is related to the input angle. Therefore the 
linear band covers the region [—2,2]. b n corresponds to 
the closed channel, whose frequency is determined by the 
interaction. The resonance occurs when the open chan- 
nel amplitude at the impurity site is zero and the local 
mode is fully excited. Therefore, we decouple (4) [tak- 
ing ao = 0] and look for localized solutions of the form 
b n = boy\ n ~ nt \ with \y\ < 1. We find that the frequency 

of this mode is uj c 



plus (minus) sign holds for e < g 2 / /3 (> g 2 //3), and 
A > (< 0) if g < e (> e). We notice that a reso- 
nance is only possible when the frequency of the open 
channel matches the frequency of the closed one. It can 
be proved that the necessary condition for having uo cc 
inside the linear band is (3 > e > 0. In any other case 
uo cc will be out of the band and no resonances will be 
observed. This condition also implies that the linear im- 
purity is absolutely necessary for having Fano resonances 
in the present model. Fig. 1(b) shows the existence region 
of the local nonlinear mode as a function of Uq. The fre- 
quency uo cc lies inside the linear band approximately for 
C/o G [1, 1.9] and, therefore, only in such a region a plane 
wave can excite a local mode and be totally reflected. 
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Fig. 2. Transmission coefficient T as a function of beam 
angle k and amplitude Uq. Bright (dark) regions denote 
high (low) T. Dashed line marks places where uok — ^cc 
and T = 0. 

The scattering problem is studied by considering an 
incoming plane wave and a localized local mode: 



b n 



% y 



g e ik(n—rii) 
p e ik(n—m) 

\n— rii I 



D e -iHn-ni) . n < n . 

; n > rtj 



By inserting this ansatz in (3) and (4) at sites n = 
Tii^rii ± 1 we find F = B + D and uok = 2cos&. 
By solving the algebraic problem we get the transmis- 
sion coefficient T = ^/B] 2 in closed form as T = 
4 sin 2 fc/[4sin 2 k + Q 2 (k)} where, 
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2A± a/4+ (e-g 2 /(3), where the 



P ±2 v /(A-cosA:) 2 - 1 

(5) 

where the plus (minus) sign is used if A — cos k > 1 (< 1). 
By fixing j3 and e, the remaining free parameters are the 
beam angle k and the amplitude Uq. Figure 2 shows the 
transmission coefficient in terms of k and Uq. First, we 
clearly see that T = (Fano resonance) appears only in 
the range Uq ~ {1, 1.9} where the frequency of the open 
channel can coincide with the frequency of the closed 
channel. Second, resonances are only possible when the 
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NLM is unstaggered (Uo > 1). This picture suggests a 
high degree of controllability of the transmission coef- 
ficient: By choosing different input powers and/or dif- 
ferent input beam angles, we can efficiently control the 
amount of light going through or reflected back from the 
impurity region. 




Fig. 3. T vs k for Uq = 0.5 (squares), 1.36 (diamonds), 
and 2 (circles). Lines correspond to theoretical T and 
symbols to numerical simulations results. 

We perform numerical simulations of (1) in order to 
confirm our theoretical findings. As an initial condition, 
we take a localized nonlinear impurity mode ?7o#' n_ni ' 
and an incoming plane wave initially centered at no: 
<t>n = <t>o exp [— a(n — no) 2 ] exp [ik(n — no)]. The initial 
amplitude </>o = 0.01 was chosen very small compared 
to L^o to fulfill the analytical criteria, no <C rii in order 
to avoid a possible initial overlap, a = 0.001 provides a 
wide spacial distribution (~ 120 sites) to correctly sim- 
ulate the scattering of a "plane wave" with a well de- 
fined k. In Fig. 3 we show results for three different values 
of [To- The agreement between the theoretical T (lines) 
and direct numerical simulations (symbols) is excellent. 
This result validates our theory and provides a good sup- 
port for observing this phenomenon in real experimental 
setups. Below the critical amplitude L^o = 1, no reso- 
nance exist because no local mode can be excited. The 
transmission profile is similar to the one found for the 
scattering of a plane wave against a nonlinear impurity 
(Uo = 0) [11]. For Uo > 1, the observation of Fano reso- 
nances is allowed. However, the stability of the nonlinear 
mode should also be considered. A small-amplitude plane 
wave can be viewed as a linear perturbation, therefore 
the localized solution has to be stable in order to numer- 
ically (and experimentally) observe the resonance. For 
Uo = 1.36, the red theoretical curve matches perfectly 
the numerical one (diamonds) where the unstaggered 
solution is stable. Stability is an extra condition that 
should be fulfilled (for nonlinear cubic impurities [7] , the 
nonlinear localized mode is always stable). Finally, for 
large amplitudes no resonances are possible because uo cc 
lies outside the band. The NLM grows in power and be- 
comes an effective wall for the plane wave. Therefore, 
total absence of transmission is expected for large Uq. 
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Fig. 4. Linear output profile for different k. The dashed 
line marks the position of the impurity site({7o = 1.36). 

Finally, for a fixed Uo = 1.36 we construct the output 
profile for different angles [see Fig. 4]. For small k- values, 
the transmission is very high achieving its maximum 
(T = 1) around tt/8. As k is increased, the transmission 
and reflection coefficients becomes of the same order of 
magnitude. For k > 7r/4, most of the energy is reflected, 
becoming a maximum exactly at the Fano resonance, 
where no transmitted light is observed. 

In conclusion, we have proposed a new possible setup 
for observing Fano resonances in optical waveguide ar- 
rays. We found a very good agreement between theory 
and numerical simulations showing that, in principle, 
this phenomenon could be observed in current exper- 
iments. We also showed a switch of nonlinear modes 
which can be controlled by changing the input power. 
We believe that our findings constitute a good example 
of the potential usefulness in using waveguide arrays for 
all-optical communication systems. 
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